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The influence of the environment in the thermal equilibrium properties of a bipartite continuous 
variable quantum system is studied. The problem is treated within a system-plus-reservoir approach. 
The considered model reproduces the conventional Brownian motion when the two particles are far 
apart and induces an effective interaction between them, depending on the choice of the spectral 
function of the bath. The coupling between the system and the environment guarantees the trans- 
lational invariance of the system in the absence of an external potential. The entanglement between 
the particles is measured by the logarithmic negativity, which is shown to monotonically decrease 
with the increase of the temperature. A range of finite temperatures is found in which entanglement 
is still induced by the reservoir. 
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I. INTRODUCTION 

The classical Brownian motion has become one of the 
most important paradigms of irreversible processes in 
physics, where out-of-equilibrium fluctuations play an 
important role. The equation that describes such an open 
system is the famous Langcvin equation, in which dissi- 
pation and noise are present. 

However, when an open system is treated quantum me- 
chanically, the problem has to be analyzed much more 
carefully. Since quantum mechanics is designed for deal- 
ing with hamiltonian systems, one is bounded to employ 
the so-called system-plus-reservoir approach to tackle 
this problem. Basically, one can consider the system and 
the environment together as an isolated system and quan- 
tize it. It is also possible, depending on the conditions 
to be obeyed [TJ EJ [3] , modelling the interaction between 
the reservoir and the system in an appropriate way and 
consider the environment as a bath of non-interacting 
harmonic oscillators. Once this is done it is possible to 
reproduce the Langevin equation for the variable which 
describes the dynamics of the system of interest in the 
adequate limits 14.. In this context of open quantum 
systems one can describe dissipation [T] and the loss of 
quantum mechanical coherence which is known in the 
literature as decoherence [5]. 

When the interest is in the dynamics of two Brown- 
ian particles in a common reservoir, some changes have 
to be made in the traditional bath-of-oscillators model. 
Those changes are described in Rcf. [6 , where the classi- 
cal equations of motion for two Brownian particles have 
been obtained. The same generalization of the model 
has been used to evaluate the quantum dynamics of two 
Brownian particles [7J. It is now our task to study the 
thermal equilibrium properties of this system. 

At first sight, one could expect that only dissipation 
and decoherence should appear in this problem. But, un- 
der certain conditions, the reservoir can correlate these 
two particles, an effect similar to the Cooper pairing in 



BCS superconductors [8]. In particular, in the model of 
Ref. [6], an effective potential between the particles ap- 
pears as a consequence of the hypothesis of the model. 
A possible correlation effect induced by this potential is 
the entanglement, a property of quantum states of great 
interest to the field of quantum information [9l [10l [11] . 
We are going to quantify the entanglement between the 
particles induced by the reservoir as a function of the 
temperature of the composite system. Some attempts to 
address this question have already been made [TH [T5] . 
but we believe that the present model [13 UJ is more ad- 
equate to account for the correct quantum behavior of 
two particles in the common heat bath. 

In section II we introduce the generalized model in 
some detail. The classical equations of motion are also 
reviewed. Then, we evaluate the equilibrium reduced 
density matrix for two quantum Brownian particles. In 
section III, we use this density matrix to understand how 
the entanglement between the particles is related to the 
temperature of the environment. 



II. MODEL AND EQUILIBRIUM DENSITY 
OPERATOR 

As we mentioned before, the possibility of modelling 
a Brownian particle within the system-plus-reservoir ap- 
proach is very well-known. In particular, if the envi- 
ronment is represented by a bath of non-interacting har- 
monic oscillators whose coordinates are coupled bilinarly 
to the coordinate of the particle of interest one can repro- 
duce the Langevin equation in the classical limit [1, 4J. 

However, depending on the specific coupling between 
the system of interest and the reservoir, a generalization 
of that conventional model is needed for the case of two 
brownian particles [B]. In what follows we review the 
important aspects of such a model for two particles in a 
common reservoir and present our main results for their 
equilibrium density matrix. This can be compared to the 
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dynamical case [7] which has been recently studied. 

We follow the Feynman path integral formalism, which 
implies that the action of the reservoir on the system of 
interest can be described in terms of the so-called influ- 
ence functional [14 . In the equilibrium case it is actually 
an Euclidean version thereof to which for simplicity we 
still refer as the influence functional. 

As an intermediate step, we show that the influence 
functional of two particles can be separated in three 
terms: one for each independent subsystem and another 
one for the effective interaction between them. It is en- 
lightening to notice that indeed the influence functional 
for each subsystem within the present method reproduces 
the result obtained by the bilinear (conventional) cou- 
pling model. 



A. Model and classical equations of motion 

The Lagrangian of the composite system reads 



L — Ls + Lj + Lr. 



(1) 



The term Ls is the Lagrangian for the system of inter- 
est, L] is the one for the interaction between the system 
and the reservoir and Lr is the Lagrangian for the bath. 
Explicitly, we have for two free particles 



L s = l -Mx\ + X -Mx\. 



(2) 



The reservoir will be described by a symmetrized col- 
lection of independent harmonic modes, 



1 N 

Lr = - TOfe \RkR-k — Lo\RkR- 



(3) 



k=l 



A non-linear coupling in the coordinates of the system 
appears in the interaction Lagragian whose form is 

Li = -~ {P-fcOi) + C- k (x 2 )]R k + [C k ( Xl ) + C k (x 3 )]R- k } ■ 
2 k 

(4) 

The local interaction of a particle with a spatially ho- 
mogeneous environment can be represented by a function 
of the type 



C k (x) = K k e 



ik'l 



(5) 



Eq. (|5| appears, for example, in the polaron problem 
[15] . This function is capable of preserving the trans- 
lational invariance of the system in the absence of an 
external potential. 



The properties of this model as well as its suitability 
to describe an effective Brownian motion are deeply dis- 
cussed in [6] and [7]. In particular, reference [6] shows 
the classical equations of motion for 



c = 



X\ + x 2 



and £ = X\ — x 2 , 



(0) 



which are respectively the center of mass and relative 
coordinates of the two particles. Those equations are 



M£(t) + (r,-rj[m)i(t) 



KffW)) 



F 6 (t) (7) 



and 



Mt{t) + ( V + rj[C{t)])t = *'dt), 



with 



mm = v 



i ±kie(t) 



and 



irkl 



Ke{t) + i 



(8) 



(9) 



(10) 



We define r\ ~ ^ fc fc 2 KfcK_fe/(/c) as the dissipation con- 
stant. The other constants are f2, a cutoff frequency for 
the oscillators of the bath, and fco, the inverse of its char- 
acteristic length scale. In fermionic environments, for 
example, this length is related to the Fermi wave number 
kF [DUE]- The functions F^(t) and F^(t) are fluctuating 
forces depending on the initial conditions [6 imposed to 
the composite system. 

Note that (|9| and ( 10 1 introduce some new phenom- 
ena into this problem. The former indicates how the net 
dissipation depends on the distance between the particles 
whereas the latter presents an effective potential between 
the particles mediated by the environment. In a regime 
in which the particles are close enough (£ << fc,^ 1 ), the 
center of mass still behaves like a Brownian particle but 
the relative coordinate describes an undamped harmonic 
oscillator of frequency ujq, with to 1 = 4r)tt/Mw. When 
the two particles are far apart (£ >> A:^ 1 ) both vari- 
ables describe standard Brownian motion with constant 
damping rj. 



B. Two particle equilibrium density operator 

The total lagrangian in Eq. ([lj generates a hamilto- 
nian H for the composite system that can be used to 
describe the density operador of the system of interest 
Ps, 



Ps = Tr R [p tot ], 



(11) 
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where 



for which the formal solution is 



Ptot 



(12) 



such that Z is a normalization constant (the partition 
function of the full system) and f3 satisfies j3 = (fcsT) -1 
where ks is the Boltzman constant and T, the tempera- 
ture, as usual. 

In the coordinate representation (2j [14] , we get 



ps(xi,yr,x 2 ,y 2> f3) = / T>[qi(r)] T>[q 2 (r)] 



<'M> ! -^{S^[qi(T)]+S^[q 2 (T)]}j x F[ Ql (r), q 2 (r)]. 

(13) 

The term Sq is the Euclidean action of the isolated sys- 
tem, 



■S;t'k(rl]= /^M^] -r ( ,-l |,/r. (14) 



with V(x) being an external potential, which is null in 
our case. The influence functional is 



P = Texp W] + ^fe(r)]} drj , 

(21) 

where Hj[g 4 (r)] = e pHR H I [q i (T)]e~ pHR and T is the 
temporal ordering operator. In what follows we explain 
the reason for using the interaction picture. 

A reservoir must be sufficiently robust so that any in- 
teraction with the system of interest is only capable of 
weakly perturbing it. Assuming such a condition, we ap- 
ply an imaginary time dependent perturbation theory for 



the bath [17]. Therefore we can expand (21 1 to second 



order in perturbation strength, which results in 



where 



1 

n 



I Y(r)dT+^ dr dT'Z(r,T'), 
Jo fl Jo Jo 



and 



y(r) = {# z [ gi (T)] +J ffzfe(r)]} 



Z{ry)^Y{r)Y{r'). 



(22) 



(23) 



(24) 



F[q 1 (r),q 2 (r)]=Ti I 



exp 



hp 



H IR {r)dT 



in which H m = H R + Hj[qi(T)] + Hi[q 2 (r)], with 



(15) 



HiMr)] = lYl {C-k[qi(r)]R k + C k [ qi {r)]R_ k } , 



(16) 



Hi[q 2 (r)] = i C -kMr)}R k + C k [q 2 (r)]R. k } 



and 



lr 1 / PkP-k . rn k oj k 



2 \ 2m k 



R k R 



k-n—k 



(17) 



(18) 



We can now identify the term e ® Hr in Eq.(19l with 



the equilibrium density operator of the reservoir, p R — 
e ~/3H a ^ an( j uge |_ ne no t a |;i on (q^ — r Yx R \p R Q ] for a gen- 
eral operator O. 



In this new notation, we see from ( 19 I that we have to 



evaluate (p). By (22 1, (23) and (24 1 we note that this 



procedure involves the evaluation of terms such (Y(t)} 
and (Z(t,t')}. Observing the definition (16) we find 



(H/[g<(r)]) =\^{C- k [q i (T)](R k ) + C k [q i (T)](R- k )}, 

k 

(25) 

which is identically zero, since (R k ) = Vfe when the 
bath is in thermal equilibrium. 



This fact implies that the functional in (19 1 can be 



equivalcntly obtained from the expansion of the following 
exponential 



It is a simple task to show that the influence functional 
can be equally written in the interaction picture as 



F[g 1 (T) I g 2 (r)]=Tt Ji [e-^«p]. (19) 
The operator p defined above satisfies the relation 



F[qi{T),q 2 {T)] = exp 



i r h P r I 2 



where 



^ = (H J [g i (r)]H J fe(r / )]). 



(26) 



(27) 



^ = -{H I [q 1 (/3)]+H I [q 2 W]}p, (20) 



As in ( 25 1 , the "tilde" has already been dropped because 



of the cyclic property of the trace. 
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Assuming again (16 1 and (17 1, we have 



A « = lT l C -k[9iir)]C- k >[q J (7 J )](R k (T)R v (7 J )) 



k.k' 



+ C- k [q l (T)}C kl [q J (T , )}(R k (T)R^ k ,(r')) 
+ C k [q l {T)]C- k \q 1 {T l )]{R_ k {T)R k ,{T 1 )) 
+ C k [q t {T)]C k . [ qj (r')](R- k (r)R-k> {/)). 



(28) 



With the help of the condition of invariance of the La- 
grangian under translation of the system of interest, only 
the terms with k — —k' survive in the sum. Then, we 
have 

k 

+ C k [q l {r)]C^ k [q 3 (r')]){R k (T)R_ k (T 1 )). (29) 

The treatment of the correlation function $fc(r — 
t') = (i?fc(T)i?_fe(r')) requires the use of the fluctuation- 
dissipation theorem |15j . which states that 



(R k (t)R- k (s)) 



(30) 



In the above relation we use X k {w) f°r the imaginary part 
of the Fourier transform of the linear response function of 
the k th oscillator. This function can be modelled in terms 
of the macroscopic behavior of the bath, since describing 
its microscopic details is not our intention here. Such an 
approach allows us to consider the harmonic oscillators of 
the environment as very weakly damped ones [B] , yielding 



X' k » = 



IkU 



m k [(uj 2 - uj 2 ) 2 + u^ 2 }' 



(31) 



which permits a simple separation of the length and time 
scales of the system. The damping constant of each bath 
oscillator is given by j k - In the low frequency limit we 
have 



X , »=/(fc)we(f2-w), 



(32) 



where we have introduced the same high frequency cut- 
off f] of the part A of this section. The Markovian dy- 
namics is achieved when we take the limit f2 — > oo and 
the function f(k) represents the nonlocal influence of the 
bath. A functional dependence like Eq.(32l for the dy- 



namical susceptibility of the bath has also been employed 

in [UGH EES]. 



In order to calculate Q k (r — r'), we make the usual 
analytical extension t = — ir and s = —ir' in (30). As 



we can see from (31 1 the response function is odd in fre- 
quency, Xfc^) = — Xfe( — which allows us to rewrite 
pOl) as 



^k(T-r') = 



cosh 










sinh 







(33) 

With the correlation function evaluated in terms of 
the phenomenological parameters of the problem, we can 
return to the influence functional and evaluate it with 
the choice we made in ((51), which leads us to 



F[qi(T),q 2 (T)] = exp 



— 2 / dr dr' a(r, t') 

h Jo Jo 



where we define 
a(r,r') = 

and 

S fe [r,r' 



k 



K k K- k [S k {r, r')$fc(r - t')\ 



2 

E 



cos[fc(g 2 (r) -Qj(r'))]- 



(34) 
(35) 

(36) 



In the beginning of section II we mentioned that we 
were going to show that the influence functional could be 
separated in three terms. Moreover, two of them could 
be interpreted as the functionals of the subsystems only. 
Now we make this statement more clear. Observe that 



(34 1 can be written in the product form F[qi(r), (^(t)] 
F[q 1 (T)]F[q 2 (T)]F mt [q 1 (T),q 2 (T)], where we define 



F[Q%(t)] = exp 



dr I dr 1 a^r, t') 
'a 



(37) 



and 



a i( T , T ') = «k«-k$k(T - r') cos[fc(g i (r) - %(>'))], 

k 

(38) 

for i = l,2. 

We can assume that the trajectory fluctuations of each 
particle are restricted to a small region compared to the 
characteristic length of the bath, fc^ 1 . Then, we expand 



the cossinc function to second order in Eq. ( 38 ) . The ze- 
roth order term will be left for the normalization of the 
density matrix. To proceed, we define 



dm u>- 



cosh 




-r'-f); 




sinh 





(39) 
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and identify, with the help of (32 1 



and 



$ fc (r-T y ) = f(k)-K(T-r'). 

Using the fact that r\ — J^k K k K -kk 2 f(k), we have 
Fh(r)} = 



(40) 



exp 



1 

~2h 



tip 



dT I dT'K{T-T'){ qi {T)- qi {T')) 2 

Iq 



(41) 

Eq.(41 1 is exactly the influence functional of one particle 



in thermal equilibrium with a bath of harmonic oscilla- 
tors with bilinear coupling [15j. 

The third term of that product is Fi n t[q±(r), qzir)], 
which involves the correlations between the coordinates 
of the two particles. However, such a correlation re- 
lates the coordinates in different imaginary times, in the 
form cos[/c(gi(r) — ^(t'))]- To deal with this problem we 
choose to apply the inverse coordinate transformation of 
that made in ([6]). That is, we write q\ and q^ in terms 
of the center of mass C and the relative coordinate £ and 
substitute them in Fi nt [qi(r), q^r)]. With some algebra 
we get 



F int [C(r),e(r)] =exp 
where we define 



-I rh,p PT 

— dr dr' A mt (T,r')K(T-r'] 

m Jo Jo 

(42) 



A mt (r, r') = K*«-*/(fc)A&!(r, r')j , (43) 

with 

A«(r,r') ee 2cos[fc(C(r) - C(r'))} cos[fc 6(r,r')} (44) 
and 



(Ty)s (m±m 



(45) 



We again evoke the fact that the fluctuations of the 
center of mass trajectory are small compared to the 
length scale of the bath, that is, C(t) p» C(t'). So, we 
expand Eq.(44l to second order and get 



A mt (r,r')= P{r,r') + Q(t,t'), 
where we define P(t,t') and Q(t,t') as 

P(t, t') = 2J2 «**-*/(&) cos[fc 6(t, t')] 



(46) 



(47) 



Q(r,r') 



k 



K k ^ k f(k)k 2 cos[k6(T,T')}. (48) 



The evaluation of the sums can be done in the same 
way as in [BJ. Briefly explaining, one can turn the 
sum into an integral by introducing the function g(k) = 
■^-KkK-kf(k) and, following [6 , we can model this func- 
tion by g(k) = e~ k / k °, where fco has been already defined. 
Surprisingly the integrals can be exactly solved and the 
result is 



and 



P(r,r / ) = --V eff (6(r,r')) 



Q(T,r') = -fj[e(r,T% 



(49) 



(50) 



with 77[0(t,t')] and V e// (6»(r, r')) defined in g and ([ToJ, 
respectively. At this point the physical meaning of the 
parameters of the model emerges. 

If we rewrite also F[qi(r)] and F[q2(r)} in terms of the 
new variables, we see that the equilibrium density matrix 
of the two brownian particles is 



psfo, = J' v[C(t)]J f V[Z(T)]exp 



t Stiff 



where 



and 



rh/3 

S eff EE / dT 
10 



mc 2 + -e 



(51) 
(52) 



1 rhP rr 

Tee-/ dr dr' v [9(r,T')}K(r-r')(ar)-ar')) 2 + 
V Jo Jo 



dr / dT'K(T-T r )(C(T)-C(r')) 2 + 
o Jo 



7T r H0 r 

rjO J dT J ^'K(T-T')V e f f [e(Ty)]- 



Tui 



dT dT'K(T-T')(Z(T)-Z(T')) 2 . (53) 



It is possible to interpret each term in the equations 
above in order to understand the solution's consistency 
and solve the Feynman path integrals. 

The center of mass £ behaves like a free quantum brow- 
nian particle in thermal equilibrium with a reservoir with 
dissipation constant ry and an "anomalous" dissipation 
given by t}[6(t, t')]. 
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The relative coordinate £ behaves like a quantum 
Brownian particle in a potential V e f / (£) in thermal equi- 
librium with a reservoir of dissipation constant r\. 

Although we have a result that agrees with the qual- 
itative interpretation induced by the classical equations 
of motion, the Feynman path integrals still have to be 
solved. An exact solution to those integrals is only avail- 
able under certain approximations such as, for exam- 
ple, the Gaussian approximation. That means that only 
terms of order 0, 1 and 2 in £, £, £ and £ in the exponential 
can be analytically computed. 

Recognizing the necessity of making approximations in 
(51 1, we will base ourselves upon physical reasoning to 



proceed further with our calculations. 

Firstly we can suppose that the interaction potential 
will be able to attract the particles. While being at- 
tracted, they experience the effect of friction until they 
end up very close to each other. So, in thermal equilib- 
rium we can assume that the reasonable regime is that 
of short relative distances and the effective potential can 
be considered harmonic with frequency ujq, which has al- 
ready been defined as u> 2 = Ar/^l/Mir. 

Besides the potential V e ff[8(T, r')], the "anomalous" 
dissipation fj[0(T, t')] also has a functional form that in- 
cludes terms of order higher than 2. We can again assume 
the regime |fco£( T )| << 1 and obtain an ohmic dissipa- 
tion, ?i[0(t,t')] ~ r/. 

These arguments transforms the formal solution in 
Eq.(51 1 into an approximate solution that carries all the 



fundamental characteristics and effects of the first one. 
Now the variables are decoupled and we can write 

Ps(Q,Q^f,^(3) = pdCf,Q,f3) p& f ,Zi,0), (54) 
where 



pdC f ,<Zi,P) = J 'z>[C(r)]exp 



and 



p 4 (£/,&,/?) = J X>[£(r)]exp 
In the equations above, we have 



1 q m 



1 <j[2] 

'n 



o[l] - 



and 

cjP] - 



M( 2 + 2 dr'K(T-r')(ar)-ar')f 
Jo 



(55) 



(56) 



dr 
(57) 



rh/3 

Jo 



4 Jo 



+ ~ote + jl <//Mr-r').i(r)-i.r'),- 



dr. 



(58) 



The solution of the thermal equilibrium density matrix 
for a Brownian particle in a harmonic potential can be 
found in [II [15] . This allows us to write 



/0c(C/,Ci,/3) = C ( exp 
and 

p${€f>€i,P) = Q ex P 



8<<Z C 2 } 



(C/ + G) 2 



2h 2 



<pf> 
2h 2 



(C/-0) 2 

(59) 



(£/-£) 2 

(60) 



The parameters of the equations are 



IjTT Jo \ (w| - uj 2 ) 2 + 4 7 j-w 



2^2 j_ 4^,2,, ,2 / V 2 



and 



(p 



2 2mjh'jj 



o \{uj 2 -uo 2 ) 2 +^ 2 oj 2 



(61) 



coth dui, 

V 2 J 

(62) 



such that 7j = rjj/(2rrij), Cj is a normalization constant 
and j — C, £■ The constants are given by 



and 



u> 2 — 0, — 2M and rj^ = 47j 



2 4^0 M . r\ 

u>c = , mt = — and ne = — . 

5 Mtt 4 2 '* 2 



(63) 



(64) 



Finally, we can return to the original variables of the 
problem, that is, the coordinates of the individual parti- 
cles. Then. 



Ps(xi,x 2 ;yi,y2,l3) = (%i E2|As|yi 2/2) 
= C c C\exp[-f(x 1 ,x 2 ,yi,y2)}, 



(65) 



where 



/(an, 2:2,2/1,2/2) 



= (2/1 ~ 2/2 + xi - z 2 ) ; 



^(P 2 ), , x2^ (2/1+ 2/2 + ^1 +X 2 ) 2 

+ 7T2"(yi ~ ^2 - ari + x 2 ) H 



2r 



32(g c 2 ) 



(P 2 



+ ^ry(yi + 2/2 -xi- X2)' 
on 



(66) 



So, we have (651 as the thermal equilibrium density ma- 



trix for two Brownian particles in a common heat bath. 
In the next section we evaluate a measure of entangle- 
ment in such state as a function of the temperature of 
the system and of the coupling between the particles and 
the reservoir. 
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III. THERMAL ENTANGLEMENT BETWEEN 
TWO QUANTUM BROWNIAN PARTICLES 

At this stage we are able to answer the question 
whether it is possible for the bath to create or keep the en- 
tanglement between two open quantum systems in ther- 
mal equilibrium with it. We describe the entanglement 
as a function of the temperature of the system and the 
coupling between the particles and the reservoir. The 
latter manifests itself through the dependence of the en- 
tanglement on the damping constant. 

We then conclude that in fact the effective potential 
mediated by the bath is responsible for the state non- 
separability. Furthermore, we see that there is a range of 
finite temperatures for which this property persists, no 
matter what the initial states of the two particles is. 

A general Gaussian state p can have its entanglement 
suitably quantified by the Logarithmic Negativity, Ejj[p] 
[9] . We see by ( 65 1 that indeed it is the case of the state 
As- 

We now make a brief review of the main aspects of this 
measure. First of all, it is worth to say that it is based 
on a criterion of separability proposed by R. Simon [TU] . 
who generalized the ideas of A. Peres [11] about positive 
partial transposition violation. 

A two mode Gaussian state can be written in terms of 
a characteristic function, defined by 



W(X) = exp 



— X a X 1 
2 



where X = {qi,Pi,q%,P2) and a is the covariance matrix 
[9]. The characteristic function relates to the density 
matrix through the relation 




0.2 



0.4 



O.d 



0.8 



1.0 



Figure 1: Logarithm Negativity as a function of the temper- 
ature A. The coupling between the system and the heat bath 
is parametrized by r. The graphic is valid for the low tem- 
perature and underdamping limit. Observe that the origin of 
the graphic is (0.1, 1.0). 



and 



2v 



-2 



2 -4detcr. 



(71) 



The new parameter z>_ is the so called sympletic eigen- 
value of the partial transposed density matrix and is re- 
lated to the logarithmic negativity [5] through 

Eh[Ps] = max[0, -lnt?_]. (72) 

After some algebra, we show that for the density ma- 



(67) trix in Eq.(65l we have 



(? f 2 > 



2h 

which implies that the entanglement is given by 



(73) 



W{X) = 



Ps 



d 2 r. 



(68) 



Eqs.(67) and (68 1 connect the state of the system with 



its covariance matrix and that is basically the reason for 
their introduction. The interest in the covariance ma- 
trix is a consequence of the fact that it caries all the 
information about the separability of the state. Its fi- 
nite dimension is fundamental for the applicability of the 
referred criteria [TUJ [11] . 

The general form of a two mode covariance matrix is 



a 7 

1 7 T 



(69) 



where a, (3 and 7 are 2x2 matrices. With them, we 
define the local sympletic invariants 



A = det a + det (3 - 2 det 7 



(70) 



■In 



'(gfKgf) 
An 2 



(74) 



To plot this function with respect to the temperature 



of the system we define, with the help of (64 1, the dimen- 
sionless variables 



r = r K 



and A = 



fk B T 



(75) 



LOq \ 

The parameter r measures the coupling between the par- 
ticle and the reservoir. The parameter A measures the 
temperature. In the low temperature (0 < A < 1) and 
underdamped (0 < r < 1) limits, the entanglement be- 
haves as in Fig.(l). 

As expected, the entanglement between the particles 
monotonically decreases with the increase of tempera- 
ture. The other significant point is that the higher the 
coupling between the system and the reservoir, the faster 
the entanglement is lost with the increase of temperature. 



8 



IV. CONCLUSIONS 

The focus of this work was the study of the thermal 
equilibrium quantum properties of two Brownian parti- 
cles in a common reservoir. More specifically, we evalu- 
ated the equilibrium density matrix for the particles and 
the entanglement between them. 

The presented model was based on the system-plus- 
reservoir approach. The excitations of the bath were de- 
scribed by harmonic oscillators and the coupling was non- 
linear in the system coordinates. The separation of the 
length scale from the time scale in the response function 
allowed us to analytically study an effective interaction 
between the particles, which is responsible for the en- 
tanglement in their state. The chosen response function 
reproduces the Brownian motion for each particle when 
they are far apart. 

We obtained a density matrix that, in terms of Feyn- 
man path integrals, presents, besides the effective poten- 
tial, the "anomalous dissipation" explicitly in its func- 
tional form. The center of mass behaved like a free quan- 
tum Brownian particle in thermal equilibrium, with a 
friction dependent on the distance between the particles. 
The relative coordinate behaved like a quantum Brow- 
nian particle in a binding potential, in thermal equilib- 
rium. 

Finally, we concluded that the environment was not 
only responsible for the loss of coherence, as intuitively 
expected, but also for the induction of entanglement. 
Even for finite temperatures this phenomenon still per- 
sists, independent of the initial conditions, since only 
thermal equilibrium was considered. The entanglement 
decreased monotonically with the increase of tempera- 
ture. Such a decrease was enforced by the increase of the 
coupling between the system and the environment. 

In our view the establishment of entanglement between 
those two particles mediated by the bath reflects the exis- 
tence of a complex ground state of the composite system 
which in general is far from being separable. 

Our results may be relevant for a better understanding 
of quantum information in non-isolated systems. They 
can also be applied to other configurations of the system 
of interest, for example, the one in which each particle is 
in an independent potential, again in the same bath. This 



model may be also adaptable to other realistic systems, 
opening the possibility to study quantum properties of 
general continuous variable open systems. 
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